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Abstract 

o . 

We study the Navier-Stokes and Euler equations of incompressible hydrodynamics 
O I in two and three spatial dimensions and show how the constraint of incompressiblility 

D ■ leads to equations of Monge- Ampere type for the stream function, when the Laplacian 

00 ■ of the pressure is known. In two dimensions a Kahler geometry is described, which is 

Q\ • associated with the Monge-Ampere problem. This Kahler structure is then generalised 

to 'two- and- a- half dimensional' flows, of which Burgers' vortex is one example. In three 
l_H ■ dimensions, we show how a generalized Calabi-Yau structure emerges in a special case. 

00 : 

1 Equations for an Incompressible Fluid 

Flow visualization methods, allied to large-scale computations of the three-dimensional 
■ incompressible Navier-Stokes equations, vividly illustrate the fact that vorticity has a 

tendency to accumulate on 'thin sets' whose morphology is characterized by quasi one- 
dimensional tubes or filaments and quasi two-dimensional sheets. This description is in 
<^ ■ itself approximate as these thin structures undergo dramatic morphological changes in 

time and space. The topology is highly complicated; sheets tend to roll-up into tube- 
^ | like structures, while tubes tangle and knot like spaghetti boiling in a pan (Vincent & 

Meneguzzi 1994). Moreover, vortex tubes usually have short lifetimes, vanishing at one 
place and reforming at another. The behaviour of Navier-Stokes flows diverge in behaviour 
from Euler flows once viscosity has taken effect in reconnection processes. Nevertheless, 
the creation and early /intermediate evolution of their vortical sets appear to be similar. 
No adequate mathematical theory has been forthcoming explaining why thin sets tend to 
be favoured. The purpose of this paper is to investigate this enduring question in the light 
of the recent advances made in in the geometry of Kahler and other complex manifolds. 
While many difficult questions remain to be solved and explored, we believe that sufficient 
evidence exists that suggests that three-dimensional turbulent vortical dynamics may be 
governed by geometric principles. 

The incompressible Navier-Stokes equations, in two or three dimensions, are 

^+u-Vu+-VP = uV 2 u, (1) 
ot p 

^ + V-(pu)=0. (2) 
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Here, u(x,t) is the fluid velocity, the pressure and density of the fluid are denoted by 
P(x,t) and p(x,t) respectively, V is the gradient operator and v is the viscosity; in the 
inviscid case when v = we have the Euler equations. The constraint imposed by the 
incompressibility condition 

V • u = 0, (3) 

is very severe. It means that the convective derivative of the density vanishes. In turn 
this means that an initially homogeneous (constant density) fluid remains constant for all 
time; p{x, 0) = p(x,t) = constant. Hereafter this density is taken as unity. Moreover, 
when (jHJ) is applied across Q it demands that velocity derivatives and the pressure are 
related by a Poisson equation 

-V 2 P = u id u^i , (4) 

where V 2 is the Laplace operator. Solving (0J) in tandem with (^Q) is computationally 
demanding and is one of the major limiting factors in how far numerical calculations can 
be driven. 

A considerable literature exists on the dynamics of vortex tubes, particularly on the 
topic of the Burgers' vortex (Burgers 1948). In an influential paper that contains substan- 
tial references, Moffatt, et al. (1994) coined the simile Burgers' vortices are the sinews of 
turbulence and thus identified the heart of the problem; that is, these filament-like vortices 
stitch together the large-scale anatomy of vortical dynamics. Despite the twisting, bend- 
ing and tangling they undergo, they appear to be the preferred states of Navier-Stokes 
turbulent flows. In fact Burgers' vortices are examples of a two- and-a-half- dimensional 
flow, which can be defined by the class of velocity fields written as (Gibbon et al. (1999)) 

u(x, y, z, t) = {ui(x, y, t), u 2 (x, y, t), zj(x, y, t)} . (5) 

This flow is linear in z in the fc-direction; thus it is stretching (or compressing) in that 
direction but is linked dynamically to its cross-sectional part. The nomenclature refers to 
the fact that it is neither fully two- nor three-dimensional but lies somewhere in-between 1 
Its components must also satisfy the divergence-free condition 

ui, x (x, y, t) + u 2 , y (x, y, t) + j(x, y, t) = . (6) 

The class of velocity fields in equation first used in Ohkitani & Gibbon (2000), is 
a more general classification of Burgers-type solutions and contains the specific form of 
the Burgers vortex solutions used in Moffat et al. (1994). Included in © are the Euler 
solutions of Stuart (1987, 1991), in which u\ and u 2 are also linear in x (say) leaving the 
dependent variables to be functions of y and t. Then stretching can occur in two directions 
thereby producing sheet-like vortical solutions. 

The differences between the three-dimensional and two-dimensional Navier-Stokes equa- 
tions are fundamental because the vortex stretching term uj-Vu in the equation for vorticity 
is present in the former but absent in the latter. Nevertheless, Lundgren (1982) has shown 
that for two-and-a-half-dimensional flows of the type 

ui = -|z7(t) + tp y u 2 = ~\yi{t) -ip x m = zi(t) (7) 

1 In the case of the three-dimensional Euler equations data can become rough very quickly; our manip- 
ulations in this paper are therefore purely formal. In fact it has been shown numerically in Ohkitani & 
Gibbon (2000) and analytically in Constantin (2000) that solutions of the type in can become singular 
in a finite time, which is consistent with observations that vortex tubes have finite life-times; the singularity 
is not real in the full three-dimensional Euler sense as it has infinite energy but indicates that the flow will 
not sustain the structure © for more than a finite time. For the possibility of a real Euler singularity see 
Kerr (1983) and Kerr (2005). 
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can be mapped into solutions of the two-dimensional Navier-Stokes equations with ip(x, y, t) 
as a stream function. 

To investigate the geometric structure behind these solutions requires certain techni- 
cal tools; these are outlined in §2 of this paper. The constraint in equation (jlj) is the 
basis of our geometric arguments, and because it is true for both the Navier-Stokes and 
Euler equations, the conclusions reached in this paper are valid for both cases. It is, 
of course, to be expected that any geometric structure should be independent of viscos- 
ity. From now on when we refer to the Navier-Stokes equations it should be implicitly 
understood that the Euler equations are also included. The Kahler structure for the two- 
dimensional Navier-Stokes equations is described in §3 and then formulated for two-and-a- 
half-dimensional Navier-Stokes flows in §4. Our results show that the necessary condition 
on the pressure for a Kahler structure to exist in two spatial dimensions (with time entering 
only as a parameter) for the two-dimensional Navier-Stokes equations is V 2 P > 0. This 
constraint is highly restrictive: by no means all two-dimensional Navier-Stokes flows would 
conform to it. More promising is the equivalent condition for two-and-a-half dimensional 
solutions of type (|7|. Theorem ^ in SI shows that these two-and-a-half-dimensional solu- 
tions have an underlying Kahler structure if V 2 P has a very large negative lower bound, 
thus associating a wide set of 'thin' solutions with the Kahler property. While the ex- 
istence of a negative finite lower bound suggests some work still needs to be done, this 
result implies that preferred vortical thin sets have a connection with a Kahler geometric 
structure that deserves further study. 

The solutions considered in this paper represent the ideal cases of straight tubes or flat 
sheets; in reality, as indicated in the first paragraph of this section, these vortical objects 
constantly undergo processes of bending and tangling. Speculatively, it is possible that 
once this process is underway, solutions move from living on a Kahler manifold in two 
complex dimensions to other complex manifolds of a higher dimension, although this is a 
much more difficult mathematical problem to address. In 33 and Appendix EI following 
Banos (2002), we describe an explicit connection between Monge-Ampere operators and 
complex manifolds of Calabi-Yau type which are associated with a restricted class of three- 
dimensional Navier-Stokes flows. This explicit connection exists only for a very special 
class of flows that require the pressure to satisfy V 2 P < 0. What this limited class of 
flows physically represent is unclear at the present time. 

The work of Roubtsov &; Roulstone (1997, 2001) showed how quaternionic and hyper- 
Kahler structures emerge in models of nearly geostrophic flows in atmosphere and ocean 
dynamics. These results were based on earlier work by Mclntyre and Roulstone and was 
reviewed by them in Mclntyre &; Roulstone (2002). It has also been shown that the 
three-dimensional Euler equations has a quaternionic structure in the dependent variables 
(Gibbon 2002). The use of different sets of dependent and independent variables in geo- 
physical models of cyclones and fronts, has facilitated some remarkable simplifications of 
otherwise hopelessly difficult nonlinear problems (see, for example, Hoskin & Bretherton 
1972). Roulstone & Sewell (1997) and Mclntyre & Roulstone (2002), describe how con- 
tact and Kahler geometries provide a framework for understanding the basis of the various 
coordinate transformations that have proven so useful in this context. This present work 
has evolved from revisiting Gibbon (2002) in the light of Roubtsov & Roulstone (2001). 

2 Differential Forms and Monge— Ampere Equations 

In this section we prepare some tools that enable us to study certain partial differential 
equations arising in incompressible Navier-Stokes flows from the point-of-view of differ- 
ential geometry. An introduction to the application of some basic elements of exterior 
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calculus to the study of partial differential equations, with application to fluid dynamics, 
can be found in Mclntyre & Roulstone (2002). Here, we shall draw largely on Lychagin et 
al. (1993) and Banos (2002). 

A Monge-Ampere equation (MAE) is a second order partial differential equation, 
which, for instance in two variables, can be written as follows: 

A(j) xx + 2B<j) xy + C(j)yy + D{4> xx 4>yy ~ <g y ) + E = 0, (8) 

where A, B, C and D are smooth functions of (x, y, tfr, (f> x , <p y ). This equation is elliptic if 

AC - AB 2 -DE>0. (9) 

In dimension n, a Monge-Ampere equation is a linear combination of the minors of the 
hessian matrix 2 of (f>. We shall refer to such equations as symplectic MAEs when the 
coefficients A,B,C and D are smooth functions of (x, y, <p x , 4> y ) £ T*M?; i.e. they are 
smooth functions on the quotient bundle J^M 2 /,/ !^, where J^IR 2 denotes the manifold 
of 1-jets on M. 2 . 

2.1 Monge— Ampere operators 

Lychagin (1979) has proposed a geometric approach to these equations, using differential 
forms on the cotangent space (i.e. the phase space). The idea is to associate with a form 3 
lo £ A n (T*M. n ), where A™ denotes the space of differential re-forms on T*IR n , the Monge- 
Ampere equation A,, = 0, where A w : C°°(R n ) -» Q n (R n ) ^ C 00 (R ri ) is the differential 
operator, A w , defined by 

A^) = (d0)*w. 

We denote by (d(j))*uj the restriction of 10 to the graph of d(p (d(p '■ M n — > T*M n is the 
differential of (f>). A form u E A n (T*IR n ) is said to be effective if uj A £1 = 0, where Q is the 
canonical symplectic form on T*]R n . Then the so called Hodge-Lepage-Lychagin theorem 
tells us that this correspondence between MAEs and effective forms is one to one. For 
instance, the Monge-Ampere equation (jHJ) is associated with the effective form 

lo = Adp A dy + B(dx A dp - dy A dq) + Cdx A dq + Ddp A dq + Edx A dy, 

where (x,y,p,q) is the symplectic system of coordinates of T*1R 2 , and on the graph of 
d(f>, p = 4> x and q = 4> y . So, for example, if we pull-back the one- form dp to the base 
space, we have dp = (j) xx dx + 4> xy dy, and then dp A dq = hess(<^>)dx A dy, where we have 
also used the skew symmetry of the wedge product. 

2.2 Monge-Ampere structures 

The geometry of MAEs in n variables can be described by a pair (0,w) G A 2 (T*M n ) x 
A n (T*M n ) such that 

1. ft is symplectic; that is, nondegenerate (f2 A SI ^ 0) and closed (df2 = 0) 

2. w is effective; that is, lo A ft = 0. 

2 We denote by hess(c^) the determinant of the hessian matrix of <j>. For example, in two variables, 
hess(0) = 4> xx (f)yy — 4>xy 

3 The use of the Greek letters w and ft is common in differential geometry; these symbols should not be 
confused with the fluid vorticity vector ui. 
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Such a pair is called a Monge-Ampere structure. In four dimensions (that is n = 2), 
this geometry can be either complex or real and this distinction coincides with the usual 
distinction between elliptic and hyperbolic, respectively, for differential equations in two 
variables. Indeed, when u G A 2 (T*R 2 ) is a non-degenerate 2-form (w Aw / 0), one can 
associate with the Monge-Ampere structure (Q,uj) G A 2 (T*IR 2 ) x A 2 (T*M 2 ) the tensor I w 
defined by 

where pf(w) is the pfaffian of ui: u A u = pf(w)(f2 A O). Thus, for the effective form lo 
associated with the MAE (J5J), the pf(w) coincides with ©. This tensor is either an almost 
complex structure or an almost product structure: 

1. A w is elliptic 44> pf(w) > 44> J 2 = -Id 

2. A w is hyperbolic 44> pf(w) < 44> I 2 = Id 
and it is integrable if and only if 

:U ) = 0. (10) 



Given a pair of two-forms (S7, cj) on T*IR 2 , such that to A Q = 0, by fixing the volume form 
in terms of Q, we can define a pseudo-riemannian metric in terms of the quadratic form 

g w {X,Y) = nXn Air (vol), X,Y <E TR , (11) 

where vol is the volume form on M 2 and tt : T*M 2 i— > 1R 2 . We can now identify our MAE 
given by u; with an almost Kahler structure given by the triple (M 2 , g^, I u ) via 

u(X,Y) = g u (I u X,Y). (12) 

One can go further and in M 4 , one can show how a natural hyper-Kahler structure emerges 
by identifying points in M. with quaternions i G H. This structure was utilised by 
Roubtsov & Roulstone (1997, 2001) in their description of nearly geostrophic models of 
meteorological flows. 



3 A Kahler Structure for two-dimensional Navier- Stokes 
flows 

We shall now show Monge-Ampere equations, in two independent variables (sat (x, y) G 
M 2 ), arise in two-dimensional and two-and-a-half-dimensional incompressible Navier-Stokes 
flows. This in turn, via the relationships described in the previous section, leads us to a 
Kahler structure for these flows. 

If the flow described by (^Q) is two-dimensional, and the fluid is incompressible, we can 
represent the velocity by 

u = kx W, (13) 

where ip(x,y,t) is a stream function and k is the local unit vector in the vertical. If we 
substitute this for the velocity in (|1J), we get 

V 2 P = -2(^4 - ^yy) (14) 
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This is an equation of Monge-Ampere type (cf. (JBJ)) for ip, given V 2 P(= P xx + P yy ), and 
it is an elliptic Monge-Ampere equation (cf. (JHJ) and @ with E = V 2 P, D = —2, A = 
B = C = 0) if 

V 2 P > 0. (15) 

Following, for example, Lychagin & Roubtsov (1983), we introduce the usual notation for 
coordinates on T*R 2 , p = ijj x , q = ip y , and then we can write (|14[) on the graph of d^ via 

u 2 d = V 2 P dx A dy - 2dp A dq; A, M = 0. (16) 

We have also on the graph of dip 

O = dx A dp + dy A dg; A n = 0, (17) 

which says simply that ip xy = ip yx . Equations (fl6|) and (|T7j) define an almost complex 
structure, I W2d , on T*M 2 , given in coordinates by 



We have 



Iij 



h 

-a 
\a J 

with V 2 P = 2a 2 . This almost complex structure is integrable (cf. Q1UJI ) in the special 



case 



V 2 P = constant . (18) 

Recall that time is merely a parameter here. When P satisfies (|18|). we can introduce the 
coordinates ^,3^, and a two-form cj^y 

<Y = x - za _1 g, 3^ = y + ia~ 1 p, tj^y = dX A d^, (19) 

then (|T1)) together with (fTTj) are equivalent to 

Au xy = 0. (20) 

To summarize, the graph of ip is a complex curve in (T*IR 2 , Iu 2d )- This is a starting point 
for a Kahler description of the incompressible two-dimensional Navier-Stokes equations. 

4 A result for two- and- a- half dimensional flows 

At this point it is appropriate to work with the two-and-a-half-dimensional Burgers solu- 
tions introduced in ^in equations ©, © and (JTJ). Based on the results of the last section, 
we shall prove a more realistic result for two-and-a-half-dimensional flows in Theorem ^ 
Lundgren (1982) made a significant advance when he showed that the class of three- 
dimensional Navier-Stokes solutions (designated in §1 as two-and-a-half-dimensional) 

ui(x,y,t) = -±j(t)x + tpy- u 2 {x,y,t) = -^(t)y -ip x (21) 

u 3 (x,y,t)=zj(t)+<t>(x,y,t) (22) 
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under the limited conditions of a constant strain 7 (i) = 70, can be mapped back to the 
two-dimensional Navier-Stokes equations under a stretched co-ordinate transformation; 
see also Majda (1986), Majda & Bertozzi (2002), Saffman (1993), and Pullin & Saffman 
(1998). In l|21j). t/j = ip(x,y,t) is a two-dimensional stream function. This idea was 
extended by Gibbon et al. (1999) to a time-dependent strain field 7 = 7(i) with the 
inclusion of a scalar (f)(x,y,t) in (J22[). The class of solutions in (|21|) . which are said to be 
of Burgers-type, is generally thought to represent the observed tube-sheet class of solutions 
in Navier-Stokes turbulent flows (Moffat et al. (1994) and Vincent & Meneguzzi (1994)). 

Depending upon the sign of ^{t) the vortex represented by (|21|) either stretches in the 
z-direction and contracts in the horizontal plane, which is the classic Burgers vortex tube, 
or vice-versa, which produces a Burgers' vortex shear layer or sheet. Thus 7, which can be 
interpreted as the aggregate effect of other vortices in the flow, acts an externally imposed 
strain function or 'puppet master', and can switch a vortex between the two extremes of 
these two topologies as we discussed in §1. 

This class of solutions is connected to the results of §2 through the following theorem, 
which is the main result of this section, and of the paper 4 : 

Theorem 1 If a two- and-a-half- dimensional Burgers-type class of solutions has a Lapla- 
cian of the pressure that is bounded by V|P > — §7 2 then any associated underlying 
two-dimensional Navier-Stokes flow is of Kdhler type. 



Proof: To prove this theorem we first need two Lemmas. Firstly let u = (ui,U2,us) be 
a candidate velocity field solution of the three-dimensional Navier-Stokes equations taken 
in the form 

u 1 = u 1 (x,y,t) u 2 = u 2 (x,y,t) u 3 = zj(x,y,t) + </>(x,y,t). (23) 

with z appearing only in U3. With this velocity field the total derivative is now 

D d d d d 

m = Wt +u 'Tx + ^d-y + ^ + ^d- z (24) 

and the vorticity vector 00 must satisfy 

-^ = Suj + uV 2 uj , (25) 

where S is the strain matrix whose elements are S{j = \ (uij + Uj t i). In the following 
Lemma v(x, y, t) = (u\, u 2 ), and V(x, y, t) is a two-dimensional pressure variable which is 
related to the full pressure P in Q31JI . The material derivative is now 

V d 

——- = — — h U • V 26) 

Vt dt v ; 



Lemma 1 (see Gibbon et al. 1999) Consider the velocity field u = (v, 2:7 + (j)); then v, 
UJ3, 4> and 7 satisfy 

— + VP = uV 2 v _i = 7W3 + u \7 2 oj 3 , (27) 

4 The notation used in this section is: V is the two-dimensional gradient and V3 is the three-dimensional 
gradient. V 2 and V| are the two- and three-dimensional Laplacians respectively (to avoid confusion with 
the symbol A in §2). 
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|£ = -7* + *>VV. (28) 

^ + 1 2 + P zz {t) = ^h- (29) 

TTie velocity field v satisfies the continuity condition diw = —7 and the second partial 
z- derivative of the pressure P zz is constrained to be spatially uniform. 

Remark: While l|27|) looks like a two-dimensional Navier-Stokes flow, the continuity 
condition implies that the two-dimensional divergence diw 7^ 0; thus an element of three- 
dimensionality remains. 

Proof: The evolution of the third velocity component U3 = jz + <f> is given by 

-P z = ^ - ^ 2 u 3 = z(^ + l 2 - ^7) + + 70 " *VV) (30) 

which, on integration with respect to z, gives 

-P(x, y, z, t) = \z 2 + 7 2 - uV 2 ^j +z(^ +1 cp- uV 2 ^j - V(x, y, t). (31) 

It is in this way that V(x,y,t) is related to P(x,y,z,t). However, from the first two 
components of the Navier-Stokes equations, we know that VP must be independent of z. 
For this to be true the coefficients of z and z 2 in (|31[) must necessarily satisfy 

^± + 7 _ „ v 2 = Cl (t) ^ + 7 2 - ^V 2 7 = c 2 (t ) . (32) 

ci (t) is an acceleration of the co-ordinate frame which can be taken as zero without loss of 
generality. Equation (|31jl shows that c 2 (i) = —P zz (t) which restricts P zz to being spatially 
uniform. To find the evolution of cu 3 we consider the strain matrix S = {Sij} 

Ul, x \ (u ljy + U2,x) 2 {zix + 4>x) \ 

I (*7„ + ^) • (33) 
H z 7x+0z) §(*7i/ + 0y) 7 / 

Working out the vorticity field to from (|2*3*|) it is easily seen that (Su)s = 7CJ3. Thus lf2lfl) 
shows that uj 3 decouples from to give the equation for uj 3 in (|27|). □ 

Now let us consider the class of Burgers' velocity fields given in 1)21(1 with a stream 
function ip(x,y,t). The strain rate variable 7 is taken as a function of time only. The 
continuity condition is now automatically satisfied. The material derivative is given by 



m = di- ^ {t) { x d- x +y dy 



Ht)[x— + y—)+ J *A^,-)- (34) 



New co-ordinates can be taken (Lundgren's transformation (Lundgren (1982))) 

s(t) = exp (J jit^dA (35) 



V2„ .-, - .1/2, 



X = S ' X w = s 



y t= s(t')dt', (36) 
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which re-scale L03 and (f> into new variables 

uz(x,y,i) = s~ l u)3{x,y,t) 4>(x,y,t) = s (j)(x,y,t). (37) 

The material derivative is 

(38) 



V 


9 


vt 


~ dt 







where ip(x,y,t) = ip(x,y,t); v = hpy,—ipxj and V = id$ + J<%. The relation between 
v = (u\, U2) and v is given by 

u l = -i 1 (t)x + s 1/2 vi; u 2 = -^(t)y + s 1/2 v 2 (39) 

and the relation between the two material derivatives in combination with the respective 
Laplacians is 

R. - „v 2 = s ( —. , - uV 2 ) . (40) 
Vt \Vt J V ; 

Introducing a new pressure variable P as 

P = s -^[V-\{x 2 + y 2 )^-\ 1 2 )} (41) 
our results can be summarized in our second Lemma: 

Lemma 2 The re-scaled velocity field v satisfies the two-dimensional re-scaled Navier- 
Stokes equations (divv = 0) 

Vv ~ ~ ~ 9 

— - + VP = vV 2 v . (42) 
Vt V ; 

The vorticity Q^(x,y,t) = —V 2 ?p and the passive scalar (f>(x,y,t) satisfy 

^ = ,V 2 u, 3 , St = „VV- (43) 

Vt Vt 

Proof: From (j35|) we note the useful result that Ds/Dt = 7s. Using ()39j) we write 

^ - = -\x (7 - |7 2 ) + s*' 2 - .V 2 ^) , (44) 

jjf " = -|„ (7 - | 7 2 ) - * 3/2 (|| - -V^ 2 ) . (45) 

Next we appeal to the definition of the pressure P in (|41[) to give the velocity pressure 
relation in 1)42(1. The results for <f> and £3 follow immediately. □ 

The proof of Theorem ^ is now ready to be completed. To obtain the full three- 
dimensional Laplacian of the pressure V 2 P we use ()41(l and ()29j) and write 



V 2 3 P = | 7 2 + S 



zA7 2 «i + — — ^ - ^V 2 v 2 



3, 



dx \Vt J dy V £>£ 

2 -s 2 V 2 P. (46) 



Thus if V 2 P satisfies the condition in Theorem^t-hen the corresponding Kahler positivity 
condition (|15[) on the Laplacian for two-dimensional flow is satisfied. ■ 

Lundgren's mapping breaks down under one condition: while the strain 7(f) can take 
either sign, if it is forever negative or for long intervals, the domain t € [0, 00] maps on 
to a finite section of the t-axis. For example, if 7 = —70 = const < then s = exp(— 70*) 
and t = Jq 1 [1 — exp(— 70^)]. Hence t G [0, 00] maps onto t £ [0, 7q" ]• 
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5 Remarks on more general 3D Navier-Stokes flows 



One of the features of this paper has been the discussion between the correspondence 
between Burgers-like flows and Kahler geometry: given the ubiquity of such flows it is 
reasonable to think of them as a 'preferred' or 'attracting' class of solutions. In ^ we 
speculated on the possibility that once the flow moves away from this class of preferred 
solutions then they may move onto complex manifolds of a higher dimension. One of 
the simplest of these is the Calabi-Yau manifold. These are smooth, compact, complex 
manifolds with a Ricci-flat Kahler metric with a holomorphically trivial canonical bundle. 

To investigate the question of other manifolds exist in more general three-dimensional 
Navier-Stokes flows we exploit the fact that the incompressibility condition allows us to 
write u in terms of a vector potential A 

u = curL4 . (47) 

We might suppose that the vector potential is some function of a scalar potential ip{x, t) 

A = (F(x, t, tp(x, t)) t G(x, t, ip(x, t)),H(x, t, tp(x, t))) , (48) 

where F, G and H are any (appropriately differentiable) functions and our notation indi- 
cates that they may have an explicit dependence on x and t. If we insert such a flow, using 
(|47|) . into (@J, then we obtain, in general, a variable-coefficient Monge- Ampere equation 
for if). As we explain in the Appendix (following Banos 2002), there is a correspon- 
dence between Monge-Ampere equations and Calabi-Yau geometries in six dimensions. 
Calabi-Yau geometries have been much studied in connection with string theory and in the 
Appendix we furnish an example of how such a geometry arises in the three-dimensional 
Navier-Stokes equations, in the case of a very special choice of vector potential. 

A further variation on this theme revolves around the addition of rotation to the 
system, which has important meteorological applications. The Euler equations of motion 
now become 

dv 

— + v-X7v + f(kxv)+X7P = 0, (49) 

where if is the angular frequency. If we examine these equations in two dimensions, with 
constant rotation, then taking the divergence of (|49[) gives 

V 2 P = ~2{lg y - ^yy) + (50) 

This is an elliptic Monge-Ampere equation if V 2 P + / 2 /2 > . The associated complex 
structure is integrable when V 2 P is a constant (cf. (|1U|) ). and in this case we can introduce 
new complex coordinates 

X = ax + i(fy + 2q), Y = ay — i(fx + 2p) (51) 

with a = (2V 2 P + / 2 ) 1/2 . Once again, (|50[) together with (j!7j) are equivalent to 

w^^dlAdy, A Wjf . =o. 

If the pressure is zero, or harmonic, then <|5U|) is suggestive of a special Lagrangian struc- 
ture. A special Lagrangian structure has also been noted in the work of Roubtsov & 
Roulstone (2001), but its role in that context is obscure (see Mclntyre & Roulstone (2002) 
equation (13.27) et seq.). 
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6 Summary 



We have shown how Kahler geometry arises in the Navier-Stokes equations of incompress- 
ible hydrodynamics, via a Monge- Ampere equation associated with (fljl. Although it is 
certainly not the case that all two-dimensional flows will satisfy the condition for the Kahler 
structure to exist, the situation looks much more promising for two-and-a-half-dimensional 
flows, of which Burgers vortex is one example. In more general three-dimensional flows, 
there is a link with geometries of Calabi-Yau type, but this connection requires further 
investigation. 

Issues relating to the existence and interpretation of Kahler structures, the integrability 
conditions, and related matters involving contact and symplectic structures, were discussed 
by Mclntyre & Roulstone (2002) in connection with various Monge-Ampere equations 
arising in geophysical fluid dynamics. The semi-geostrophic equations of meteorology, 
which are a particularly useful model for studying the formation of fronts, were the starting 
point in Mclntyre & Roulstone op. cit. for an investigation into the role of novel coordinate 
systems, similar to those we have found here in (|19|) and (j51|) . In semi-geostrophic theories, 
such coordinates facilitate significant simplifications of difficult nonlinear problems, and 
they are associated with canonical Hamiltonian formulations of these systems. Issues 
relating to contact and symplectic geometry may also be relevant to the results presented 
in this paper, and this suggests one direction for further study. 

Finally, it is perhaps appropriate to make a comment here about compressible flows. 
In two dimensions, a velocity potential x( x ^) ma Y De introduced, so that in and © 
we have u(x, t) = Vx + k x Vt/>. We may take the divergence of (^Q) and this yields an 
equation for ijj of the form 

H^ly ~ ^xxlpyy) + ^Xxyi^xx ~ Ipyy) = R, (52) 

where R is a function of p and P, and their spatial derivatives; it is also a function of the 
first spatial derivatives of ip, and it will also depend on the derivatives of x with respect 
to space and time (e.g. Xx, Xxy, Xxxt and Xyyy etc.). In other words, we have a variable- 
coefficient Monge-Ampere-type equation for t/j, when p, P and x (and their derivatives) 
are known. It is therefore possible, in principle, that a Kahler structure may be associated 
with this type of equation; however, in practice, such a structure will be difficult to identify. 

A Calabi— Yau structures and 3-D-Navier-Stokes flows 

In six dimensions (T*M 3 ), there is again a correspondence between real/complex geometry 
and nondegenerate Monge-Ampere structures (Banos 2002). A Monge-Ampere structure 
(n,uj) G A 2 (iV 6 ) x A 3 (iV 6 ) on a 6-dimensional manifold N (N = T*M 3 ) is said to be 
nondegenerate if uj is nondegenerate in the sense of Hitchin (2000): the Hitchin pfaffian 
X(u) £ C°°(N) never vanishes (cf. lf55 |) -([57 jl below). We note that: 

1. A(a>) > if and only if uj is the sum of two decomposable forms, that is 

u> = ax A o<2 A a 3 + 0i A 02 A (3% 

2. X(uj) < if and only if uj is the sum of two decomposable complex forms, that is 

uj = (a\ + i0\) A (a 2 + %0i) A (a 3 + i/3 3 ) + (a\ - i0\) A (a 2 - ifh) A (a 3 - i0 3 ). 

Hitchin has associated with such a nondegenerate form a tensor K u : N — > TN <g> T*N 
which is either an almost complex structure if A(cj) < or an almost product structure if 
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A(w) > 0. Moreover, Lychagin and Roubtsov (1983) have defined a metric g w on N which 
is compatible with K u and which has signature (3, 3) if \(u>) > and (6, 0) or (4, 2) if 
X(uj) < 0. So, one can associate with a nondegenerate Monge- Ampere structure (Q,u;) on 
a 6-dimensional manifold an almost Kahler structure (which is "real" if A(w) > 0) whose 
Kahler form is £1. Moreover this almost Kahler structure is "normalized" by the two 
decomposable three forms associated with w. we use then the terminology of generalized 
Calabi-Yau structure (Banos 2002). 

For example, the Monge- Ampere structure associated with the "real" MAE in three 
variables (x, y, z) 

hess(0) = 1 
is, in the coordinates (x,y, z,p,q,r), the pair 

I Q = dx A dp + dy A dq + dz A dr 
I oj = dp A dq A dr — dx A dy A dz 



and the real Kahler structure on T*R 3 is 



The Monge-Ampere structure associated with the special Lagrangian equation 








C 3 : 



hess(0) 
is the pair 

\ Q = dx A dp + dy A dq + dz A dr 

\uj = 9((dx + idp) A {dy + idq) A (dz + idr)) 

and the underlying Kahler structure is the canonical Kahler structure on T* 



9u 



K„ 



It is important to note that the geometry associated with a MAE A w = of real type 
(A(o>) > 0) is essentially real but it is very similar to the classic Kahler geometry. In 
particular, when this geometry is integrable, there exists a potential <I> and a coordinate 
systems (xj,Pi)j=i ; 2,3 on T*R 6 such that 




and 



9u 



det 



d 2 § 

dxidp. 



d 2 $ 



■dxi ■ dpj 



f(x)g(p). 



K dxjdpj 

We have observed that the (elliptic) Monge-Ampere equation plays an intriguing role in 
incompressible two-dimensional and the two-and-a-half-dimensional Navier-Stokes equa- 
tions, and that this is made explicit via the introduction of the stream function. As we 
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remark in for more general three-dimensional Navier-Stokes flows, a stream function 
can play a role as suggested in (|48j): here we examine perhaps the simplest choice 



(53) 



where now ip = tp(x, y, z, t) (the minus signs are used simply to show the relationship with 
the two-dimensional case viz. v = k x VV> = V x (0, 0, —ip)). Note that (|53j) is simply one 
choice that leads to some interesting features as described below. This choice of vector 
potential is, of course, very special, and merely serves to illustrate the possible connection 
with generalized Calabi-Yau structures. 

If we substitute (|4"Tj) with (|53~|) for u into (J3J) (and set p = 1), and use the notation 
ip x = Pi/^y = 9; = r ; then we find that we have another Monge-Ampere equation that 
can be written as 



cu 3d = V 2 Pdx A dy A dz + 



2 dx A dp A dr + 2 dr A dg A dz + 2 dg A dy A dp 
2 dy A dq A dr - 2 dp A dg A dx - 2 dr A dz A dp 
2 dx A dq A dr - 2 dp A dy A dr - 2 dp A dt? A dz; 



A W3d = 0. 



(54) 



It is worth mentioning that since the canonical symplectic form Q = dx A dp + dy A dq + 
dz A dr vanishes on the graph of ip, one can replace co by u + 9 A f2, with any 1-form 
on T*]R 3 . Following Lychagin & Roubtsov (1983), we choose among all these forms, the 
unique one which is effective (its product with f2 is zero). This form is 

u)q = u> — -(A-0j) A U, 



with 



^A^ + ^A^+^A^ 
dx dp dy dq dz dr 



j u. 



Three invariants are associated with the effective form (|54j) : the Lychagin-Roubtsov metric 
Q, the Hitchin tensor K and the Hitchin pfaffian A (see Hitchin (2000) and Banos (2002)). 
They are defined by 



Q(U,V) 
A 



~± 2 ((U^o) A(Fj^o)) , 

^ikQkj > 

~ Tr(K 2 ) . 



(55) 
(56) 
(57) 
(58) 



When A < 0, the tensor J 



describes a Calabi-Yau geometry. We obtain here A 
with a complex geometry if and only if 



K is an almost complex structure and the pair (Q, J) 
128V 2 P: our equation is associated 



V 2 P < 0. 



(59) 



The Lychagin-Roubtsov metric is 



Q = 16 



(a 2 


a 2 


a 2 








o \ 


a 2 


a 2 


-a 2 











a 2 


-a 2 


a 2 























1 


1 











1 





-1 


\o 








1 


-1 


0/ 



(60) 
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and the Hitchin complex structure is 



J = 



1 



/ 





1 1 \ 

10-1 

1-10 

-a 2 

r 2 

a 2 / 



(61) 



V2a 




2 



,2 



a 



with V 2 P = —4a 2 . Once again, this complex structure is integrable in the special case 
when V 2 P is constant. 
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